Generalized Green Functions and current 
correlations in the TASEP 

A.M. Povolotskyi't, V.B. Priezzhev^'*, G.M. Schiitz^'H 
^ : January 18, 2011 



l> 






> 



m 

o 
o 



cd 



Abstract 



We study correlation functions of the totally asymmetric simple 
exclusion process (TASEP) in discrete time with backward sequen- 
tial update. We prove a determinantal formula for the generalized 
Green function which describes transitions between positions of parti- 
cles at different individual time moments. In particular, the general- 
ized Green function defines a probability ineasure at staircase lines on 
the space-time plane. The marginals of this measure are the TASEP 
correlation functions in the space-time region not covered by the stan- 
\^ I dard Green function approach. As an example, we calculate the cur- 

rent correlation function that is the joint probability distribution of 



On . times taken by selected particles to travel given distance. An asymp- 



totic analysis shows that current fluctuations converge to the Airy2 
process. 
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1 Introduction 

The totally asymmetric simple exclusion process (TASEP) is one of 
the most studied models of non-equilibrium statistical mechanics [STJ 
\TT\ \T7\ [30] . A rough guide in this field would distinguish between 
continuous time and discrete time processes. The latter are subdi- 
vided into several groups in accordance with updating rules: parallel, 
sequential, sublattice parallel, etc.[24J. The detailed microscopic in- 
formation about the time evolution in the TASEP is provided by the 
transition probability or the Green function of the master equation. 
The Green function describing the continuous time evolution of A^ 
particles on the infinite one-dimensional lattice has been obtained in 
the form of a determinant of a A x A matrix in [29]. For periodic 
boundary conditions the Green function is given by the A^-fold sum of 
similar determinants [20j . The determinantal expression of the Green 
function for the parallel update has been found in [22j by the Bethe 
ansatz techniques and in [4J by induction in time. The ring geometry 
for the parallel update was considered in }23j . The model we deal with 
in this paper is the discrete time TASEP with sequential update. A 
formal derivation of the determinant for the Green function for this 
process has been undertaken in [21] by a geometrical treatment of the 
Bethe ansatz and in [25J directly using the Bethe ansatz. 

The Green function defines a probability measure on the set of 
particle configurations. The correlation functions, i.e. distributions 
of positions of selected particles, are marginals of this measure. Re- 
cently, considerable progress has been achieved in the calculation of 
correlation functions. In the pioneering work by Johansson, [T^, the 
distribution of the position of a tagged particle at a given time was ob- 
tained for the TASEP with parallel update and step initial conditions. 
The technique used involved several mappings: from the TASEP to the 
last passage percolation problem, then to the problem of the longest 
nondecreasing subsequence of a generalized permutation, and, finally, 
via Robinson-Schensted-Knuth's correspondence, to the statistics of 
pairs of Young tableaux. In contrast, the same result can be obtained 
by direct summation of the Green function over the whole set of fi- 
nal configurations constrained within a suitably chosen domain. In 
this complementary way, Johansson's result was generalized to the 
TASEP with backward sequential update in [25j and extended to dif- 
ferent initial conditions for the continuous time TASEP in fT9]. The 
distributions were presented in the form of Fredholm determinants 
with polynomial kernels and were proved to converge asymptotically 
to the universal Tracy- Widom distributions of the largest eigenvalues 
in the Gaussian random matrix ensembles, unitary and orthogonal for 



step and flat initial conditions respectively. The Tracy- Widom distri- 
butions were shown to be universal scaling functions of the Kardar- 
Parisi-Zhang (KPZ) universality class [16j. 

The direct summation of the Green function hardly helps for the 
multi-particle distributions. Instead, it was observed in [27] that the 
determinantal formula for the Green function dependent on N particle 
coordinates can be treated as the marginal of an auxiliary determinan- 
tal measure on much larger space of A^(A^ — l)/2 variables, which are 
identified with coordinates of fictitious particles with free fermionic 
interaction. The marginals of this measure are calculable using the 
construction of the determinantal point process, from where the cor- 
relation functions of interest follow. Similarly to the one-particle dis- 
tributions, they can be represented in the form of Fredholm determi- 
nants and which asymptotically converge to the universal scaling limits 
defining the Airyi, Airy2, e.t.c. processes, which were claimed to be 
the hallmarks of the KPZ class. The results for joint distributions of 
positions of several particles at a given time have been obtained in 
O [21 [5] for different updates and initial conditions. The reduction to 
a suitable determinantal process turned out to be efficient also for the 
calculation of more general correlation functions [8j. In particular, the 
chain of mappings similar to that from the Johansson's work was used 
in [13] to find the auxiliary determinantal measure, a.k.a. the Schur 
process, which allowed a calculation of the distribution of positions of 
a tagged particle at different moments of time. 

The most general correlation functions, which refer to the posi- 
tions of selected particles at selected, though generally different, time 
moments, were obtained in [H [1]. The authors considered a clus- 
ter of particles initially compact, hopping to the right. Its evolution 
was sliced into as many slices as the number of selected time mo- 
ments, each slice connecting two subsequent time moments. Then the 
correlation function was constructed as a convolution of the Green 
functions, where each Green function described the evolution within 
a single time slice and fixed the position of another selected particle. 
Relying on the fact that in the TASEP the rightward motion of par- 
ticles is independent of the particles on the left of them, each Green 
function described only the evolution of particles on the right of the 
particle position fixed at the previous step. This structure stipulates 
a weak ordering of the time moments within the correlation function: 
the more to the left is the initial position of the particle, the earlier is 
the corresponding time moment. Here "weak" refers to the fact that 
the times (particles) selected for different particles (times) are also 
allowed to be equal, which reproduces the known cases of fixed time 
(tagged particle) correlations respectively. This ordering was named 



space- like [H [T] . Such a multilevel construction gave rise to even more 
complicated auxiliary determinantal point process than for the fixed 
time case, which finally resulted in the correlation functions on space- 
like point configurations. Of course, the limiting processes obtained 
from the asymptotical analysis of correlation functions at the space- 
like paths have the expected KPZ-specific form. 

The space-like ordering used in ^11] is a technical condition orig- 
inating from the method used. Therefore, one expects that it would 
not be a constraint for the KPZ universality. Indeed, basing on the 
observed slow decorrelation phenomena, it was proved in |121 [10] that 
the limiting processes are the same along any space-time path, except 
special directions given by the characteristic lines of the hydro dynamic 
conservation law associated to the large scale particle dynamics. How- 
ever, no microscopic results on correlation functions beyond the dis- 
tributions of particle positions on the space-like paths were available 
up to date. 

The aim of the present paper is to extend the set of known cor- 
relation functions by using a generalization of the conventional Green 
function, which describes transitions between space-time point con- 
figurations where different particles are associated with different time 
moments and different space positions. Such a generalization was 
proposed in [9]. It was shown by means of a trajectory treatment of 
the Bethe ansatz that, under special conditions on space-time points 
of initial and final configurations, the generalized Green function has 
a determinantal representation similar to that of the regular Green 
function. The idea now is to supply the generalized Green function 
with a meaning of probability measure on some set of space-time point 
configurations passed during the TASEP evolution, and to calculate 
the TASEP correlation functions as the marginals of this measure. We 
partially implement this program here. We define the space-time point 
sets on staircase lines named boundaries by requiring them to satisfy 
two conditions. First, the Green function describing the transition 
between point configurations within these sets is of the determinatal 
form. Second, the TASEP evolution has a Markov property at the 
boundary needed to claim that the probability of a space-time tra- 
jectory of the system factorizes into the product of probabilities of 
the two parts obtained by cutting the particle trajectories along the 
boundary. Due to this fact, the required probability measure can be 
defined as the probability of configurations of space-time points from 
where the particle trajectories go after they have traversed the bound- 
ary. Remarkably, the time coordinates of points in the configurations 
available at the boundaries are ordered in an opposite way than those 
in the space-like configurations, and, hence, we have a new class of 



configurations to deal with. Tliey, liowever, do not provide a full 
complement to the set of space- like configurations from ^ [1] , because 
there is also a restriction on the space coordinates of the points, which 
are required to be strictly decreasing with the number of a particle. 

After discussing the general construction, giving the probability 
measure on the new class of configurations, we still need to present a 
determinantal process to proceed with the correlation functions. The 
further analysis is limited to a particular example of the generalized 
boundary, where N particles exit from N fixed subsequent space points 
at arbitrary time. For the step initial conditions, this suggests that 
all particles make the same number of steps. Then, the correlation 
functions of interest are the joint distributions of times taken by se- 
lected particles to travel a given distance. For this case we obtain the 
corresponding determinantal measure and complete the calculation of 
correlation functions. We should note that similar auxiliary determi- 
nantal measure for the same problem in the case of parallel update 
was proposed in |13j, where the term current correlations was used. 
Therefore, we keep on the same terminology here. The auxiliary de- 
terminantal measure was obtained in [T3] using the above mentioned 
chain of mappings. This method is restricted to the step initial con- 
ditions only. At the same time, the generalized Green function give 
a more general tool applicable to a wider class of initial conditions 
and finial configurations. Therefore, the results of the present paper 
on current correlation functions should be viewed as a test example of 
the method used. The more general cases will be considered elsewhere. 

The paper is organized as follows. In Section 2 we define the 
discrete time TASEP with sequential update and present some basic 
properties. Then we proceed with the derivations of the Green func- 
tions. To be self-contained, we start with a regular derivation of the 
Green function for sequential update, using an inductive solution of 
the Chapman-Kolmogorov equation. This is done in Section 3. In the 
Section 4, we introduce the generalized Green function. In the first 
part of this section. Proposition 13.81 we present a rigorous proof of the 
determinantal formula first proposed in [9j. The formula is restricted 
to the case when initial and final point configurations are of the special 
form, which we call admissible, defined by the ordering of coordinates 
in time and space. In the second part of Section 4 we introduce the 
boundary sets, where the admissible configuration live. We prove that 
the generalized Green functions define the probability measure at the 
boundaries and possesses the Markov property, implying that the con- 
volution of two Green functions performed over the boundary yields 
another Green function. These results are summarized in Proposi- 
tions 14.91 In Section 5, we use the Sasamoto's formalism of auxiliary 



variables |27] to rewrite the generalized Green function in the form of 
a signed determinantal point process and write its correlation kernel 
in the Proposition I5.4[ The main results of the paper on the current 
correlations are given in Section 6. Theorem 16. II gives the joint distri- 
bution of particle currents in the form of Fredholm determinant. The 
asymptotic form of this determinant and convergence to the Airy2 
process is the content of Theorem 16.31 The latter is given without 
a detailed mathematical proof. Only the sketch of the saddle-point 
analysis of the kernel is presented in Lemma 16.21 A conclusion and a 
review of perspectives are in Section 7. In particular an application of 
the generalized Green functions to more general situations is briefly 
discussed. 



2 Discrete time TASEP 

Consider A^ particles on the one-dimensional integer lattice. A config- 
uration of the system 



takes its values in the set 



xgZ^ (2.1) 



{{xi,...,xn) G Z^ : xi > X2 > ••• > xn} , (2.2) 



which is an A^-tuple of strictly increasing coordinates of particles 
{xi, . . . ,xm)- The strictly increasing order reflects the exclusion con- 
dition that two particles cannot occupy the same site. 

The Markovian dynamics of this process for a single discrete time 
step can be defined as follows. Starting from the rightmost particle at 
position xi < oo, we let each particle jump to its right neighbouring 
site with probability p, provided that the target site is empty. If 
the target site is occupied, the probability for a particle to stay is 
1. Therefore the TASEP is a random process which is given as a 
sequence of configurations x^,x^, . . . ,x^. We refer to such a sequence 
as a trajectory of the system up to time t. Every trajectory is realized 
with probability 

P{x°, ...,x*) = Pi{x^\x^-^) . . . Pi{x'^\x^)Pi{x^\x'^)Po{x'^), (2.3) 

where Po{x) is the initial probability of configuration x. The one- 
step transition probability Pi{x\y) from configuration y to x takes 

the form 

N 

Pi{x\y) = Y{0{xi-yi,Xi-i-yi), (2.4) 

i=l 



where 

d{k,l) = {q + pSi^i) Sk,o + p6k,i, (2.5) 

and we formally define xq = oo. The parameter p, the hopping prob- 
ability, varies in the range 

0<p<l, (2.6) 

and we also define 

q = l-p. (2.7) 

In other words, we exclude the trivial deterministic limiting cases 
p = 0, 1 from our consideration. The above transition probabilities 
correspond to the discrete time dynamics with backward sequential 
update, when at each time step the position of the first particle is 
updated first, then of the second particle, etc. Since the hopping is 
totally asymmetric, the first A^ particles do not "feel" the presence 
of any particle to their left. Hence, when considering the dynamics 
of the first A^ particles, one may study without loss of generality the 
system with exactly A^ particles initially. Then for any fixed t the 
number of configurations and also the number of possible trajectories 
(with non-zero probability) is finite. 

3 Green function 

Consider the probability for the system to be in a configuration x after 
t time steps, t > 0, 

Gt{x\y)= Y. Pix'',...,x'), (3.1) 

{a;0,a;i,...,a;*-i} 

given 

Po(3J°)=5y,,o. (3.2) 

The sum is taken over all trajectories, the final configuration being 
fixed 

a;* = X. (3.3) 

The transition probability satisfies the recurrence relation 

Gt{x\y)= Y, Gi{x\x')Gt-i{x'\y), (3.4) 

which is a particular case of the Chapman-Kolmogorov equation. It 
follows directly from the definitions (|3.1|) and (|2.4|) that 

Go{x\x°) = 6^^^o (3.5) 

Gi{x\x') = Pi{x\x'). (3.6) 



An explicit form of Gt{x,y) for arbitrary time t can be found as a 
Green function (GF) of the equation (|3.4p that satisfies ()3.5p . (|3.6p . 
The solution is particularly simple in the case of the TASEP, having 
the form of a determinant of an A^ x A^ matrix. A first determinantal 
formula for the GF Gt{x,y) has been found by Schiitz [29] for the 
continuous time TASEP. A similar formula for the TASEP with back- 
ward sequential update has been obtained by the use of trajectory 
treatment of the Bethe Ansatz |2Ij. As pointed out in [22] it can be 
proved rigorously using the determinantal approach of [29]. For the 
sake of self-containedness and as introduction to the approach used 
further below for multi-time correlations we present here a complete 
proof of this formula by induction. 

Proposition 3.1 The GF of J^'.-^p has the determinantal form 

Gt{x\y) = det[Fj_i{xi - yj, t)]i,j=i,...,Ar, (3.7) 

where 



"^'^'*^" 0, t<0 



(3.8) 



and the integration contour T encircles the origin, while the point w=l 
stays outside. 

Proof. The proposition can be proved in three steps. First we prove 
that the r.h.s. of (|3.7p satisfies the initial conditions (|3.5p . then that it 
satisfies the recurrence relation ()3.4p , and finally show that it correctly 
defines Gi{x\x') according to (|3.6p . so that it can be used to obtain 
Gt{x\y) for t > 1. 

We first expand the determinant representation of Gt{x\y) given 
in r.h.s. of ()3.7p into a sum over permutations: 



N 

Gt{x\y)= 5^(-l)l'^IJ]i^..-i(^.-ya,,i), (3.9) 

(J€Si^ 4 = 1 

where a = (cJi, . . . , ajy) is a permutation of N integers (1, . . . , N) and 
the summation runs over all elements of the symmetric group Sn ■ To 
show that Gt{x\y) satisfies the initial conditions we need to evaluate 
this sum in the particular case t = 0. To this end, we refer to the 
first statement of the Lemma lA.ll proved in appendix \^ It directly 
follows from this statement that the only summand which can be 
nonzero in the r.h.s. of ()3.9p in the case i = is the one consisting 



of the diagonal elements of the matrix [Fj^i{xi — yj,0)]^^^ •<^. Then 
the initial conditions are straightforwardly verified, 

Goix\y) = 5^,y. (3.10) 

This completes the first step of the proof. 

To prove the recurrence relation (|3.4|) for the r.h.s of (j3.7|) we use 
recurrence relations [^ for the functions Fn{x,t): 

Fn{x,t) = qFn{x,t- 1) +pFn{x - l,t - 1), (3.11) 

Fn{x + l,t) = Fn{x,t) - Fn^l{x,t). (3.12) 

Using ()3.1ip we can expand the determinant in the r.h.s. of ()3.7p into 
a sum of similar determinants, where we have either -Fj-j (xj — yi,t — l) 
or Fi-j{xj — iji — l,t — 1) instead of Fi-j{xj — yi,t) with coefficients 
q and p respectively. In other words, we obtain another recurrence 
relation 

Gt{x\y)= Y, G'r{x\x')Gt-i{x'\y), (3.13) 



similar to (j3.4p , but with transition probabilities for the last step writ- 
ten for N non-interacting particles: 



N 



Gf'=^(a:|^') = n^o(^^-^^'l) (3.14) 



While the configuration x is supposed to be from Z^, given by (|2.2p . 
the range of summation in ()3.14p is the set 

Z^ = {{xi,...,xn) GZ^:xi >X2 > ••• >X7v}, (3.15) 

which admits coincidence of coordinates of two successive particles, 
x[ = a;^_|_i. Of course, the quantity Gt-i{x \y) being defined for x' E 
Z^ does not have the meaning of a probability in Z> . Nevertheless 
we can extend formally the determinant ()3.7p to this set. 

Now we show that ()3.13p with ()3.14p is equivalent to (|3.4p . Com- 
pare first quantities Gi^^{x\x ) and Gi{x\x ). Let us consider a sum- 
mand of ()3.9p corresponding to a particular permutation cr at t = 1. 
The second statement of the Lemma [A. II restricts the range of permu- 
tations a as well as the range of configurations x' for this summand to 
be nonzero at given x. To characterize the structure of permutations, 
we use the notion of cycles. Given the permutation a the cycle is a 
set of elements (ii, . . . , i^) such that 

i2 = ai,... ,ik = CFi^_^,ii = ai^. (3.16) 



(a) 



\i x^e 



X3 X2 x^ 




X3 X2 



(b) 




Xo X2 X, 



Figure 1: Clusters Ck, (a) and D^, (b) for A; = 7 



Any permutation can be decomposed into a set of disjoint cycles. 
According to the Lemma lA.ll the permutations contributing to the 
sum can contain cycles of two types only. First, these are trivial cycles, 
at = i, which imply either x^ = Xi 01 x^ = Xi — 1. They contribute 
to the product the terms -Fo(0, 1) = q and Fq{1, 1) = p respectively. 
Second, the cycles which consist of the indices of successive particles 
packed into compact clusters: 



ii 



X, 



Xi 



k + l 



(3.17) 
(3.18) 



,ik = i -\-k — 1 
I 

■ 1 Xi-\-k~l = Xj^j^f^^i 

for 1 < i < A^ — 1 and 2< A;<A^ — i + 1, so that positions of clusters 
are the same in x and x' . The latter case gives the terms of type 

\k-l ■ 



Cu 



-1 



-^i2-h l-^ii 



Xi^A) 



^ik-ik-i\Xik-i 



2^ifc) t)-rij_jj.(Xjj. 



Xi,A) 



(3.19) 
which is illustrated by the diagram in FigjT^. As a result, we obtain 



Ck = {-ir~'Fi{l, lf-'F_k+x{-k + 1, 1) = p'^-'q. (3.20) 

The transition probability Gi'^{x\x') in ()3.13p also contains the 
terms Fo(0, 1) = q and Fo(l, 1) = p corresponding to free particles 
staying or making a jump with x^ = Xi 01 x^ = Xi — 1 respectively. 
In addition one must take into account cases when two particles from 
the same site of x' are coming to neighboring sites of a;. In general 
they can also belong to a bigger cluster of particles ahead making a 
step as well. This process, shown in FigjTb, contributes the term 



Dk = F^ixi - Xj, 1) . . . Fo(xi 



+k-2 



i+fc_2,l)-^o(2;j+fc-l 



(3.21) 



10 



where 



Xi-\-\ — Xi 1, . . . , Xi-|_/^_i — Xi fC -f- i, 

•^i ^ Xi — i, . . . , Xi^f,_i ^ Xi — k -\- L, [6.22) 



so that 



Dk = Foil, l)''^Fo(0, 1) = p'^-^q (3.23) 

Thus, there is a one-to-one correspondence between Gi^^{x\x ) with 



X G Z> and the nonzero summands entering the determinant expan- 



sion (I33D for Gi{x\x') with cc' G Z^. 

Consider now the terms Gt~i{x'\y) in the r.h.s. of ()3.4p and ()3.13p 
for the configurations x' , which differ by groups of arguments defined 
by p.lSp and (|3.22p respectively, i.e. 

Gt-i{... ,Xi,Xi-l,...,Xi- k + l,...\y) (3.24) 

and 

Gt-i{...,Xi -l,Xi-2,. . . ,Xi- k + l,Xi- k + 1,. . . \y). (3.25) 

The first k — 1 arguments in (j3.24p are shifted by one with respect 
to those in (|3.25p . Using the identity ()3.12p and performing column 
operations in (j3.7p for columns of the matrix Fj-i{x^ — yj\t — 1) con- 
taining shifted arguments, we equalize the determinants for ()3.24p and 
(j3.25p and therefore prove that (j3.4p and (|3.13p have the same form. 

As the final step, we show that the determinant formula for Gi(a;|a;') 
complies with the definition ()2.4p of Pi{x\x'). According to (|2.4p a 
particle which jumps a step forward brings a factor of p into Pi{x\x'). 
A particle which does not jump occupies the same site in x and in 
x' and can bring either the factor of q or the factor of 1. The lat- 
ter occurs when there is a particle in the right neighboring site that 
keeps its place as well. In other words, an isolated compact cluster 
of particles keeping their positions when the system jumps from x' to 
X always brings only a factor of q, no matter how many particles are 
in the cluster (by isolated cluster we mean that it does not belong to 
any bigger cluster). 

We need to check if the r.h.s. of ()3.9p taken with t = 1 gives the 
same result. One can see that permutations contributing nontrivially 
to the sum in (|3.9p permute only the indices of subsequent particles 
constituting compact clusters, which keep their positions unchanged. 
Therefore the sum (|3.9p can be factorized into factors each coming 
from an isolated cluster and the factors of p coming from every single 
jumping particle. The latter is already in correspondence with (|2.4p . 
What we need is to sum up all the contributions from an isolated 



11 



cluster. The result is q which can be proved by induction. Indeed, for 
a cluster consisting from a single isolated particle there is only a term 
i*b(0, 1) = q corresponding to the trivial cycle. Let us suppose that 
it is valid for all clusters consisting of i < fc particles. For k particles 
the last particle can enter the cycle of the type ()3.17|3.18p , shown in 
the FigH^, together with {i — 1), 1 < i < k, particles next to it. The 
contribution of such cycles is Cj while the sum over all contributions 
of the other (k — i) particles brings the factor of q by the induction 
hypothesis valid for i < k. An exception is i = k, when no particles 
in the cluster remain anymore, resulting simply in C^- Summing up 
all these contributions we have 

k-l 
qY,Ci + Ck = q, (3.26) 

i=l 

which completes the induction. ■ 

4 Generalized Green Function 

A different view at the TASEP evolution is a representation of the par- 
ticle trajectories as an ensemble of non-intersecting lattice paths. This 
is convenient because then we can consider a wider class of boundary 
conditions, i.e., initial and final configurations which refer to different 
times for different paths. 

4.1 Ensembles of Lattice Paths 

Consider the two-dimensional integer lattice L = T? . A point of this 
lattice is associated with two integers (x,t), x corresponding to a par- 
ticle location and t representing discrete time. Let us fix two points 
{x^,t^) G L and (x,t) G L, such that t > to- 

Definition 4.1 A directed lattice path of the length {t — to) with the 
starting point {x^,t^) and the endpoint {x,t) is a subset of L of the 
form 

i=0 

where 

x*-*° = X, (4.2) 

f = t° + i, (4.3) 

and X* G Z can take any values for i = 1, . . . ,t — t^ — 1. 
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A path is associated with a weight W ( H) A ) defined as a prod- 
uct of weights of elementary (two-point) weights, 

t-to-i 



w 






n W[{{x\f),{x'+\f+')}]. (4.4) 



1=0 



In the definition of lattice paths, the adjective "lattice" refers to the 
connectivity structure defined by a given set of nonzero elementary 
weights, which can be represented as bonds and diagonal links on a 
2D lattice with vertices in L. The paths which have nonzero weight, 
which are the only ones of interest, can be drawn on this lattice. In our 
case, the lattice paths are intended to model the stochastic trajectories 
of particles. Correspondingly, the set of nonzero elementary weights 
is defined by the allowed particle jumps. For example, for the ID 
directed Bernoulli random walk, which is a particular case of the above 
dynamics with a single particle, the nonzero elementary weights are 

W {{{x, t), (x, t + 1)}) = q, W ({(x, t), (x + 1, t + 1)}) = p, (4.5) 

where (x,t) G Z^. 

Within these definitions, the GF of a single particle for the directed 
Bernoulli random walk, which is the N = 1 case of the GF defined in 
the previous section, is given by the sum of the weights of all possible 
paths with given starting and end points: 

Gt{x\x')= Yl ^[nSifT- (4-6) 

To interpret the results of previous section for the A^-particle TASEP 
GF, we have to consider A^ interacting paths on the lattice. Note 
that the results for the GF connecting two particle configurations at 
different times strictly depend on the fact that the particles are non- 
colliding and stay ordered during the whole evolution. This is also 
true for their space-time trajectories, i.e. for the corresponding lattice 
paths. Indexed at the initial time moment, the points associated with 
different paths have the same order at any later time. Our aim is to 
define a generalization of the GF in such a way that the initial and 
final positions of different particles would refer to different time mo- 
ments. For this definition to be of practical use, we would like to keep 
the path ordering condition. It turns out, that this is still possible if 
we restrict the range of initial and final configurations: 

Definition 4.2 A set of N indexed points of L 

{X,t) = ((xi,ti),...,(xAr,tAf)) 
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is called admissible if the following conditions are satisfied simultane- 
ously 

xi > X2 > ■ ■ ■ > xn (4.7) 

and 

ti<t2<--- <tN. (4.8) 

Definition 4.3 Given two admissible N -point configurations 

(^°,t°) = ((x?,t?),...,(x^,4)) (4.9) 

(a;,t) = ((xi,ti),...,(xjv,t7v)) (4.10) 

such that 

ti>t'^i for i = l,...,N, (4.11) 

an N-path is a subset of L x {1, . . . , N} of the form 

nSif = {(ni)j::f ;,..., (n^);f-f-;}. (4.12) 

The number 

T = tN-t\ (4.13) 

is referred to as the length of the N-path. 

The interaction in tiie TASEP implies that the statistical weight 
W [n] of an AT-path is not in general the product of the weights of A*" 
single paths. Namely, it has a form of the product of terms depending 
on two adjacent paths, 

Af-l 

W[U] = W [Hi] J] W [Iii+i\Iii] . (4.14) 

i=l 

The weight W [nj+illlj] is a product of two point terms for the second 
path, rtj+i 

i^[n,+i|n,]= \{ w\{{xi^,,t^^^,{x^,t^\R, 

i=o 

(4.15) 
conditioned by the configuration of the first path IIj. All the nonzero 
two-point terms are defined as follows 

«'[{(x,*).(. + M + i)}|ni = {^; IXZX'i^^A- (""> 

Now we are in a position to introduce the generalized Green func- 
tion (GGF). 
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Definition 4.4 Given two N -point configurations {x,t) and {xo,to) 
satisfying the conditions of definition \4-3\ a GGF from (a;o,to) io 
{x,t) is 

G{{x,t)\{xo,to))= V wln^P] (4.18) 



,.O.tO^ 



(X,t) 






where summation is over all possible N -paths between {x^,t^) and 
{x,t). 

Obviously the GF Gt {x\xo) defined in the Section [3] is a particular 
case of the GGF with t^ = ■ ■ ■ = t% = t° andti = ■ ■ ■ = t^ = t+t°. Of 
course, it does not depend on t^ because of the time shift invariance. 
The main result of this section is the proof that the determinantal 
form (|3.7p of Gt {x\xo) can be extended to the GGF. 

Proposition 4.5 Given two admissible]}] N -point configurations {x,t) 
and {xo,to), such that ti > t? for i = 1,. . . ,N, the GGF can be ex- 
pressed in the determinantal form: 

G((a;,t)|(a;o,to)) = det[F_,_i(xi-x°,ti-tO)]i,j=i,...,^, (4.19) 

where Fn{x,t) is defined in i!^.8\} . 

Proof. Consider an A^-path 11; A ' . It follows from the definition 
of weights, that if we cut an A^-path along points of equal time, the 
parts obtained will be independent of each other. To be specific, let 
us fix t', ti < t' < t]\f. In general, A'^ paths forming the N- path can 
be subdivided into three groups: ending up before or at t'; starting 
before t' and ending later than t'; starting at or later than t'. In other 
words, there are two integers, 1 < A;i < A;2 < A^, such that 



h, ■ ■ ■ ,tki < i 




(4.20) 


tki+l,- ■ ■ ,tk2 > *' > *fci + l)--- ;*fc2 




(4.21) 






(4.22) 


Let us fix two pairs of configurations of k2 points 






(^o,iV = ((x?,t;),...,(xO„4)) 




(4.23) 


{X',t')^ = ((xi,tl),...,(Xfc,,tfcJ,(x',^ + i,t'),-- 


, i^'k. 


t')()4.24) 



^ As shown in [2, the range of initial and final configurations, which ensure the validity 
of the proposition, is even broader. Why we limit our consideration by the admissible 
configurations only will be clear from the subsection 4.2. 
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and (A^ — ^2) points 



(x\t') 



{x,t)^ 



((xfci+i, tfci+i), . . . , (xAT, tiv)) (4.26) 



where (x'^ _(_^, t'), . . . , (x'^ , t') are the intermediate points of the paths 
Ilfc^+i, . . . ,11^2 respectively. For the weights being nonzero they have 
to be bounded to 

Xk,+i > 4i+i > • • • > 4^ > <> (4-27) 



which ensures that the configurations (|4.23|) - (|4.26|) are admissible. 
Therefore, we can represent our A^-path as a union 



i(x°,t°)- 



ix',t'), 






(4.28) 



of fc2-pa'th and (A^ — /i;i)-path respectively, which are clearly indepen- 
dent, so that the weight of the A^-path is equal to the product of two 
weights: 



W 



n 



{x,t) 



w 



n 



(a,',t')- 



w 



n 



{x',t')+ 

(x,t)+ 



(4.29) 



Then, we can split the summation over A^-paths into independent 
summations over all possible realizations of Tl,,\,\ and 11^^ 's , 
each resulting in the corresponding GGF, and a summation over the 
intermediate points x'j^_^_^_^, . . . , x'^^ which vary in the range ()4.27p : 

G{ix,t)\{x^t°))= (4.30) 

^ G{{x,t)^\{x',t')^)G{{x',t'U{x°,t»)_) 



{<, 



■KJ 



Using this formula and Proposition 13.11 we can construct the GGF 
recursively, attaching one time step at a time. 

First step: Consider an A^-path of unit length T = tjy — t^ = 1. 
In this case all the paths constituting the A'^-path must fall within the 
same unit time interval, 



'-1 

h 



tN 



t" 

iO + 1. 



(4.31) 
(4.32) 



Obviously, a GGF in this case coincides with the ordinary GF dis- 
cussed in Section [3l 



G((a;,t)|(a;°,0)=Gi(a;|;r°), 



(4.33) 
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which fits ()4.19p . thus proving the proposition for this special case. 

Recursion: Suppose ()4.19p is valid for the A^-paths of any length 
up to (t — 1). Consider A^-paths of length t and apply (|4.30p with 
t' = t — 1, dividing an A^-path of length t into parts of lengths t — 1 
and the unit length. Again, the unit length part is given by the regular 
GF: 

G{{x,t) + \{x',t')+)= (4.34) 

Gi ((xfci+i, . . . ,xn)\{x'j._^^i,. . . ,x'j^^,Xf^^^-y, . . . ,Xj^)) , 

Then ()4.30p can be rewritten as 
G{ix,t)\ix°,t°)) = Yl G((^',0-|(^°,t°)-) 

xGi {{xk,+i, . . . ,X7v)|(4i+i> • • • >4^)) ^x;^2+i,x0^+i • • •'^x'^,x5^(4-35) 

where we replaced the starting coordinates x^ ,^, . . . , x^ by summa- 
tion variables a^fc^+i, ■ • • ,x'j^ respectively, fixing their values with the 
corresponding Kronecker delta symbols. 

Here we can formally release the summation domain ()4.27p to the 
whole Z> ~ ^ , as one of the factors of the expression under the sum 
is zero beyond this domain. The first factor together with Kronecker 
deltas can be packed into a single determinant in the domain ()4.27p . 

= det [F.^M - x^,t[ - t°)] f,^/4.36) 



where we set 



x^ = xi,...,x,^^=Xk^, (4.37) 

ti = ti,... ,tj.^ = tk^,tj^^_^^ = • • • = tjv = t , (4.38) 

This can be verified by noting that for 

l<i<k2, k2 + l<j<N, (4.39) 

the following matrix elements are zero 

Fj_M-^",ti-t^j) = 0. (4.40) 

Indeed, we have either 

t'i - t° < (4.41) 
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or 



t[-t°j = 0, (4.42) 

X- - x° = x'i- x'j > 0, (4.43) 

as the configuration (x',t')^ is admissible. Then (j4.40p follows from 
the properties of Fn{x,t), (|A.6IIA.8"]l . Thus the matrix on the r.h.s. of 
()4.36p has a block form 



where 



B 



A = [Fj.,{x', - x^t', - t])]^.^^^^,^^ (4.45) 

= \F,-U-^A-m^,=,,^^,...,N (4-46) 



By definition 

detA = G((a;',0-l(^°,*°)-) (4.47) 

and from Proposition 13.11 



5„, „o ... (5^/ ^0 . (4.48) 



detB = Go((xfc2+i,...,x'^)|(x^2+i'---'^A^)) 

Then ()4.36p follows from 

det(^ ^j=detAdetB. (4.49) 

It remains to prove that a substitution of the determinant for 
G\ (x\x') and the r.h.s of (|4.36p into (j4.35p results in the determinant 
(|4.19p . This is completely analogous to the proof of the recurrence re- 
lation in Proposition 13.11 up to the fact that only the matrix elements 
from the columns with the indices j = /ci + 1 , . . . , A^ are involved in 
the summation. ■ 

4.2 Probabilistic content 

The GF introduced in Section [3] defines a probability measure on the 
set of particle configurations (|2.2p and thus has the natural interpreta- 
tion of a conditional probability to be in configuration x at time step 
t given an initial configuration x^ . However, the GF is also a marginal 
measure defined on the set of all TASEP realizations with exactly t 
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steps. When we sum over all possible realizations provided the final 
configuration at step t is fixed, we get the Green function Gt{x\x^). 
Moreover, due to the Markov character of the process, it obviously 
has the same meaning also on the set of process realizations of any 
longer duration oi t' > t steps. Then, the GF still gives a probability 
of a particle configuration to appear at the intermediate step t. Inte- 
gration of the GF over a subset of configurations, where the positions 
of some particles are constrained to a given domain, yields a corre- 
sponding new marginal probability, also called correlation function. 
An objective of the present article is an extension of this scheme to 
marginal probabilities referring to particle positions at different times. 
To this end, we aim to treat the GGF as a marginal measure defined 
over more general sets than the set of fixed time particle configura- 
tions. Before discussing the general scheme, we first outline it with 
the example of a single particle. 

4.2.1 Single particle (iV=l) 

Consider a set of lattice points 

B = {ixi,ti)},<,z, (4.50) 

where for any i £ Z either 

(xi+i, tj+i) = {xi + 1, ti) (4.51) 

or 

(xi+i, ti+i) = {xi,ti - 1). (4.52) 

We refer to ^B as a boundary set implying that it can be viewed as a 
boundary of the set 

L= [J {{x',t') eL:x' <x,t' <t,}. (4.53) 

{x,t)eB 

Note that 

BcL. (4.54) 

The boundary set has the following obvious property 

Property I. Any directed Benoulli random walk starting inside L 
and ending in its complement, L'^ = L\L, crosses the bound- 
ary. Furthermore, once having left L it can never return. More 

specifically, consider a path HL A such that [x ,t ) G L and 
(x, t) £ L'^. Then, the path necessarily contains exactly one point 
{x',t') G n[f;P, such that 

{x',t')eB (4.55) 
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and 

{x",t")eL^ (4.56) 

for any {x",t") £ njf;*°^ with t" > t' . 
For each {x, t) £ B we define an exit probability 

Pf((x,t))= Yl W[{{x,t),{x',t + 1)}] 

= Yl W[{{x,t),{x',t')}]- (4.57) 

{x',t')eL'= 

A generahzation of the probabihstic meaning of the GF follows 
from the following simple claim. 

Proposition 4.6 Given {x^,t^) G L, the quantity 

g«((x',t')l(x°,t°)) = P,^iix',t'))Gi{x',t')\ix^,t^)) (4.58) 

defines a probability measure on B. 

Proof. Suppose first that L is bounded from above in time direction, 
i.e. there exists a time t such that 

LcLt = {{x', t') £L:t' < t, x' G Z}. (4.59) 

Then any path with starting point (x*^, ^) G L and end point at (x, t), 
with any x G Z, definitely leaves L via a unique point of B due to 
property I. 

Consider a set of paths of length T = t — to 

nf''°^ = {ng'J) : X G Z}, (4.60) 

with all possible endpoints. Note that the weights we consider main- 
tain the probability conservation 

Y H^[n] = j;G((x,i)|(xo,to)) = l. (4.61) 

Therefore weights H^[n/i^'V ] define a probability measure on Q,^ ' 



assigning probability to each particular path. We can now decompose 

(•3,0 ^0-1 

any path from il^ ' into three parts 
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where {x',t') £ B and {x",t' + 1) G L^. The Markov property imphes 
that the weight of a path is a product of the weights of its parts : 

Using this property we can sum the weight of ()4.62p over all paths 
from ilj ' with {x',t') £ B fixed and {x",t' + 1) varying in L^: 

g«((x',t')l(x°,t°))= E WiieqllM) 

!.n(£n[''°'*°^:{x',t')&B;{x",t'+i)eL'=j 

(4.64) 
The summation over the last segment yields 1 due to ()4.6ip . while the 
rest results in the r.h.s. of (j4.58p . 

Therefore, Q {{x',t')\{x^,t^)) gives a marginal probability for a 

path form Q^' ' to exit from L at a particular point {x',t') of the 

)(^'o,*o) 



"t 
boundary. As any path from QY''^''''^' definitely crosses the boundary, 
this in particular implies 



Y, g^{{x',t')\{x',t'))= Yl ^[n] = l. (4.65) 

ix',t')€B {nenl"0''o^} 

Now one can release the condition (j4.59p which ensures that every 
path from Oj^°' ° crosses the boundary before time t. This may not 
happen if the boundary approaches infinite time at some finite value 
of the coordinate x > x^, which is to say that 

BD[J{x,t* + i) (4.66) 

for some t*. Then, for any t there exist paths in Oj^"' " which do 
not reach the boundary within (t — to) time steps. A simple estimate, 
however, shows that their total measure in Q^' ° vanishes as i — )• oo. 

Prob (n C 0f°'*°) : np|i3 = 0) = O (p^-^^g*-*") . (4.67) 

Therefore we have 

Y: e^((x',oi(x°,t°))=iim Y g^{{x',mx',t')) 

(x',t')GB {x',t')eBf]Lt 



lim V W[U] = 1(4.68) 

|ncn["O''o^nns5^0} 



t—^00 



This proves the proposition. ■ 
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The restriction of the set of paths to ilj^"' ' ahows us to inter- 
pret the function Q {{x,t)\{x^,t^)) in terms of the particle evolution. 
Specifically, Q,^' ° can be treated as the set of all possible sequences 
of (t — t^) positions of a particle that started its walk at site x^. Cor- 
respondingly, Q {{x,t)\{x^,t^)) gives a marginal probability over this 
set for the sequence, in which space-time position of the particle leaves 
the domain L right after taking the value {x, t) € B. We may conclude 
that when {x,t + 1) ^ L, Q^{{x,t)\{x^ ,t^)) gives the probability for a 
particle to be in site x at time t, irrespectivly of what the next step 
is. When (x, t + 1) G L, this is a probability for the particle to jump 
out of site X at time t. 

The decomposition ()4.62p and the Markov property (j4.63p allow 
us to obtain a formula 

G{{x,t)\{x^,t^))= Y. G{ix,t)\ix",t' + l)) 

{x',t')eB;{x",t'+l)eL'' 

xW{{{x',t'),{x",t' + l)})G{{x',t')\{x^t'')). (4.69) 

for the convolution of the GGF associated with a general boundary 
which generalizes the formula (j4.30p we used for convolution along 
fixed time configurations. 

We consider two examples of boundaries and corresponding Green 
functions we deal with below: 

Example 4.7 Fixed time boundary, 

^(•'*) ={(2;,t) :xeZ} (4.70) 

with t G Z being fixed. This is a case already discussed in section ^. 
Its Green function is given by 

g^'-"((x,t)|(x°,t°)) = G((x,t)|(x°,t°)) =p— %*-*"-+^° (^*:^°o 

(4.71) 
ifO< (x-x°) < (t-t°) andG{{x,t)\{x^,t'^)) = otherwise. 

Example 4.8 Fixed space boundary, 

i3(^'-) = {(x,t):tGZ} (4.72) 

withx > xq & Z being fixed. At each point the GF Q^ ' {{x,t)\{x^ ,t^)) 
has to be supplied with a "jump off the boundary" probability p: 

g^'"'((x,t)|(xO,iO))=pG((x,t)|(xO,tO)), (4.73) 
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with G{{x,t)\{x^,t^)) like in ( [-<^.7Jp . Then the normalization and con- 
volution read: 

f; jf-^^^^q^-'^-^^^^ (^ - '^\ = 1, (4.74) 

t=x-xO+tO ^ ^ 



E (:::i:)(:::'r\)^(::::).(^-^'.*-. 



t—x+x' 

E 

tl=t0+x'-X0 



4.2.2 A^ particles 

We generalize now the concept of a boundary set to the A^-paths. 
To motivate our approach we first observe the following. A natural 
candidate for an A^-boundary would be an arbitrary collection of A^ 
boundaries of the form (|4.50p . 

B = {Bi,...,Bn} (4.76) 

However, an attempt to do such a straightforward generalization fails. 
There are two obstacles. First we note that the factorization of the 
weight of an A^-path into the product of weights of its parts does not 
take place for arbitrary endpoints of these parts. Indeed, consider two 
adjacent steps of particles i and i — 1 belonging to the A^-path 11 

n D {(x, t), {x, t + 1)}, \J{{x + 1, t), {x + l,t + !)},_! (4.77) 

Within n the weight associated with {(x,t), (x,t + l)}j is 

W ({(x, t), (x,t + l)}|{(x + l,t), (x + l,t + 1)}) = 1. (4.78) 

On the other hand if we separate a part 11 of 11, 11 C 11, so that 

n D {(x, t), (x, t + i)}„n ;^ {(x + 1, t), (x + i,t + i)}i_i, (4.79) 

a similar weight within 11 will be 

W{{{x,t),{x,t + l)h) = q. (4.80) 

Thus, if we use the decomposition separating 11 from 11, the weight 
is not factorized and hence we cannot use a convolution formula like 
()4.30p . In general, the difficulty occurs when the endpoints of segments 
{(x, t), (x, t + l)}j and {(x + 1, t), (x + 1, t + l)}i-i occur in the sets 
Li and L'^_i respectively. While they interact in the common A^-path, 
their weights become different when they do not "see each other" in 
different subsets. 
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The second difficulty is tliat the GGF we dealt with is defined 
only for the case of N- paths with admissible configurations at the 
endpoints. At the same time, mutual positions of the endpoints of 
different paths within an A^-path varying in arbitrary boundaries can 
violate obviously the admissibility constraint. 

We need to find such a form of an A-boundary that it would ensure 

(a) the weight factorizability when the A^-path is divided along B and 

(b) admissibility of initial and final configurations for any A^-path with 
nonzero weight. Both requirements can be fulfilled with A-boundaries 
constructed as follows. We introduce a translation operator T^ trans- 
lating a set of space-time points by n steps in x direction. 

Tn\J{{x^,ti)} = [J{{xi + n,ti)} (4.81) 

i i 

Now take an arbitrary boundary B of the form (j4.50p . The A- 
boundary we are looking for can be constructed from ;B as a collection 
of A^ copies of B each translated one step back with respect to the 
previous one: 

B={B,T^iB,...,T^N+iB) (4.82) 

Indeed, this construction excludes the possibility to have the endpoints 
of segments {(x, i), (x, t + l)}j and {{x + 1, t), {x + l,t + l)}j-i in the 
sets Li and L1_^ simultaneously. In addition, as soon as path IIj leaves 
Bi at {xi,ti), the coordinates of point (xj+i,ij+i) of Bi+i where path 
IIj+i leaves i3j+i obey ij+i > ti and x^+i < Xi due to the noncrossing 
conditions ()4.16p . (|4.17p . The latter inequalities are the admissibility 
conditions. 

Using these properties we can straightforwardly generalize the one- 
particle scheme to A particles, replacing all the arguments above to its 
bold font counterparts. As a result, we have the following statement 
which provides us with a probabilistic interpretation of the GGF. 



Proposition 4.9 Let B he given by \4.82 ). A quantity 



g^{{x,t)\{x',t')) = Pf{{x,t))G{{x,t)\{x\t% (4.83) 

defines a probability measure on the set of admissible configurations 
on B, and the GGF satisfies the following recurrent relation. 

G{{x,t)\{x\t')) = Y. G{x,t\x",t") 

{x',t')eB;{x",t")eL'= 

X W{{{x', t'), {x", t")})G{x',t'\x°, t°)(4.84) 
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5 Signed determinantal process 

The notion of GGF allows us to consider the TASEP realizations, 
where different particles enter and end up their evolution at different 
moments of time, such that their initial and final configurations are 
admissible and the latter belong to an A^-boundary of the form (j4.82p . 
The aim of this section is to show that for specific initial conditions the 
measure over the realizations of these processes defined by the GGF 
according to Proposition 14.91 can be represented as the conditional 
measure of certain signed determinantal point process, which allows 
for the calculation of the multi-point correlation functions. Below we 
give one example of this calculation for the TASEP evolution of A^ 
particles starting at the same moment of time t = at sites x^ = 
l-i,i = 1,2,. ..,N. 

Consider a set of the A^-paths starting at ((0, 0), . . . , (1 — A^, 0)) and 
terminating by making a step away from A^-boundary B constructed 
as A^ copies of the fixed x boundary (j4.72p : 

N 

B = [j{{x-i + N,ti):tieZ}, (5.1) 

j=i 

where the subscript i refers to a boundary set associated with a path 
IIj within the A^-path corresponding to i-th particle. We are interested 
in the joint distribution 

P = Prob ({t„, < ai} f]{tn, < as} fl • • • f|{i„„ < a„}) (5.2) 

of times t„^ , . . . , t„^ , when m < N particles labeled by indices rii < 
77,2 < • • • < iT-m jump off positious Xni , • • • x^^ taken from an array Xi = 
x + N — i,i = 1,. . . ,N and tm < in2 ^ " " " ^ ^n^- Similar calculations 
for other cases of a wider class of initial and final configurations, which 
also can be performed by the GGF technique, will be done elsewhere. 
Consider an auxiliary signed point process over the subsets of 
Z>a; X {1, . . . ,n} of the form 

T= \J K,<T^_i,<...,<rncZ>(,_i)X{l,...,n} (5.3) 

l<n<Af 

given by the measure 

•^('^) = ^ n det[c/>„(rf,r7+i)]^+i,det [^f (t^-JI^Iq , (5-4) 

^ n=0 

where we define the functions 

U^.y) = {l:lf, (5.5) 
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and 

*f (i) = {-lfF^k{x + N-k-l,t), (5.6) 

where 

F„(x,t) = -^/ ^(g + ^)*(l-u;)-V, (5.7) 

The integral representation holds for t G Z. This is unlike Fn{x,t), 
which coincides with Fn{x,t), when t >0, and vanishes at t < 0, see 

We also introduce fictitious variables t"~^, 1 < n < A^, which are 
effectively less than any rj^, so that 

'/'n(r„\i,r;+i)=p (5.8) 

for J = 1, .., n+1. If we consider r" as coordinates of fictitious particles 
at the n-th time step, then T^_^_l corresponds to a particle entering into 
the system from a reservoir on the left |27j . 

Having defined the measure Ai {T) we are able to interpret the 
GGF in the following useful way. 

Proposition 5.1 Given 

{x^t^) = ((0,0),...,(-iV + l,0)) (5.9) 

{x,t) = {{x + N-l,ti),---,{x,tN)) (5.10) 

with ti < t2 < ■ ■ ■ < t]\f £ Z>a; and x > —N + 1 the GGF associated 
with the boundary Ii5.1\) is a marginal of the measure M 

(5.11) 

To prove this statement one represents the GGF as a sum over 
the auxiliary time variables in a way similar to that used for space 
variables in j27[ HI El |3] and for time variables in |19j . To this end we 
first state and prove the following lemma. 

Lemma 5.2 Given initial and final configurations 115. 9\) and 115. 10\) 
respectively, the GGF G{{x,t)\{x^ ,t^)) can be represented as a sum: 

G{{x,t)\{x^,t^)) = {-P)^ 

x^ det[F_Af+i+i(x + i,r^i)]ij=o,...,Af_i (5.12) 

D 

where t( = tj, j = 1, . . . ,N, and the summation variables take their 
values in the domain 

D = {ri£ Z>.,2 < . < J < iV|r/ > r/-\r/ > rZ+^j. (5.13) 
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Proof. Let x{a, b) be the summation variable with its upper and 
lower limits a, b. Then the order of summation in ()5.12p is given by 
the sequence 

N( N-l_-. N N (^N-l_. ^N-U N~l( N-2 _-. s 
IN\'N-1 ^^-^J^ ' N-1\' N~2 ^^ ' N-1J^ ' N-1\' N-2 ■'-i-^A---) 

^N f^N-l _ 1 ^N-l\ N^l, N~2 _ -, ^N-2\ N-2fN-3 _ i ^\ 
^N-2VN-3 ^^'N-2)^'N-2\'N-?, ^^'N-2)^'N-2VN-3 -"-i-^j)---) 
N( N-l_-, ^N-U fc/ fc-l_i ^ 

r2^(rf-i-l,rf-^),...,r|(rii-l,x). 

The variables {r", 1 < i < n < N} and inequalities between them are 
shown in Fig. [2] for the case A^ = 4. The total number of summations 
is iV(iV - l)/2. 

The summation formula for the functions Fn{x,t) is 

p^Fn{x,t) = Fn+lix + l,t2 + l)-Fn+l{x + l,ti) (5.14) 

t=ti 

The first summation with t^ converts the last column in the de- 
terminant from F^N+i+iix + i,T^) into F_N+2+i{x + 1 + i,rj!^zl), 
i = 0,1, . . . ,N — 1. Note that generally the summation ()5.14p pro- 
duces two terms. However the second one, F_7v+2+i(2; + 1 + i,x), al- 
ways vanishes. The second summation with tJ^_i produces two terms 

F^N+2+iix + 1 + i, rl^Z2) - F-N+2+i{x + l + i, t^zI) in the previous 
column. The second term coincides with the last column and disap- 
pears. The third summation with t^Zi converts the last column into 
F_Ar+3+j(x + 2 + ZjT^Jg)- Continuing, we notice that each row of 

summations with r^* for fixed k increases the index n in Fn{x,t) by 
one in each column starting from (N + 1 — A;)-th one and increases 
argument x by one in these columns. Performing all summations, we 
obtain the determinant 



det 



F^N+i{x, rf ) ... Fo{x + N- 1,tI) 
Fo(x + iV-l,rf) ... Fn-i{x + 2N-2,tI] 



(5.15) 



As the variables r^ are nonnegative by assumption for j = 1, . . . ,N, 
we can omit the tildes, to obtain the determinant coinciding with 
G{{x,t)\{x^ ,t^)) up to a factor (— l)^^/^! after reversing of the order 
of columns. ■ 

Now we are in a position to prove the previous proposition. 
Proof. Proof of Proposition [XT] As is illustrated in Figl21 the inequal- 
ities in ()5.3p follow directly from the inequalities in ()5.13p . Then, the 
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X 



X+1 



X+2 



X+3 



Figure 2: The set of auxiliary variables for A^ = 4. Open circles are inte- 
gration variables, closed circles are fixed moments of time t^ = t^, r^ 



r? 



4, 'i — ^3, 

^2, tI = ti. Arrows show inequalities between variables. Red lines 
correspond to variables with a fixed upper index and can be considered as 
discrete-time steps of an auxiliary process. 
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statement of the proposition follows from Lemma [5. 2 1 and the fact that 
^ J.1 r N{N+i)/2 n-aD 

The latter can be expressed in the language of symmetric functions 
(see [I8] for notations and the facts used). Specifically, up to the 
coefficient of p the function (pnix,y) is nothing but the complete sym- 
metric function specialized at single parameter equal to 1, (j)n{x,y) = 
phy^xi^)- Then, by the generalized Jacobi- Trudy formula 

det[</>„(rf,T;+i)]i<,,,<„+i=p("+i)s^,„+i)/;,(„)(l), (5.17) 

where we defined partitions 

Xin) =^n^l_x^^n^2-X...,T^+n-X,0,... 

ior n = 1, . . . ,N. Here s ^^(n) / x(n-i) (x) is the skew Schur function. For 
X = 1, it is equal 1 provided that A'"-* D A'""^-* and the skew partition 
(a'"' — A^"^^-*) is a horizontal strip. If at least one of these conditions 



is violated, s^{n)/^(n-i){l) = [T5]. Translated back to the variables 

r/, these two conditions exactly coincide with the inequalities ()5.13p 
defining the domain D. ■ 

The measure Ai (T) has the same functional form as in [271 [HI [3] . 
In particular, Lemma 3.4 of [3j can be directly applied here. There 
is a difference in the form of functions ^^ (t) and in the space where 
the variables r- live, which is Z^x rather than Z. This difference does 
not affect the applicability of the Lemma, which is formulated in a 
rather abstract fashion, though has to be taken into account when 
obtaining the final expressions. According to Lemma 3.4 the multi- 
point correlation functions of A^(7~) are determinantal. To write down 
the kernel we introduce 



{(t>ni * <t>ni+l * ■■■ *(t>n2-l){x,y), "-i < 722 

(5.18) 






where (a * b){x, y) = Y.zeL^^ a(^' ^)^(^' V)^ and 

*^,(T) = (<^"'^*^j^_,.)(r). (5.19) 



Consider functions 



{(00 * </'(''"))(r?,r), . . . , (0„_2 * 0--^'")(r-r2,r),,/.„_i(rr\r)}. 

(5.20) 
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They are linearly independent and hence can serve as a basis of an 
n-dimensional linear space Vn- We construct another basis of Vn, 
{<I>^(r), j = 0, . . . , n — 1}, which is fixed by the orthogonality relations 

Y, ^nr)^]ir) = 6,,. (5.21) 



Then, under the 

Assumption (A) : 0n(T^+i,''") = c„$o(t) with some c^ / 0, n 

1 N, 



L, . . . , X T , 



the kernel has the form 

n2 

K{n,,n;n2,T2) = -<P^-^''^-\n,r2) + J2K[-kin)KU{r2). (5.22) 

k=l 
Let us now make this construction explicit. 

Lemma 5.3 The functions ^^(r), *J??(t) have the following integral 
representation. 

ii,n(^r) = i ^(l-p^^^^\\yj-l)kuj-+N-k-2 (5 23) 

/ro.i 2vr? V w J 

^n^r) = pi ^(l-p^-Zl\ " (t; - 1)-^- V-^--ti24) 
■' hi 27rz V V J 

The contours of integration Fq,!, Ti encircle the poles w = 0, 1 and v = 
1 anticlockwise respectively, leaving all the other singularities outside. 

Proof. To construct ^^(t) we start with \I'^(r), /c > 0, which is 
given by ()5.6p . and follow recursively the definition (|5.19p 

^1-\W=pE^^(*)- (5.25) 

t=T 

Note that for /c > the expression under the integral ()5.24p for ^^(r) 
has no pole at w = 1, and therefore the integral complies with the def- 
inition ()3.8p of the function F„(x, r) and the definition ()5.6p of ^^(t). 
Then, substituting (j5.24p into (j5.25p . we can interchange integration 
and summation. To this end, we must ensure that the resulting sum 
converges uniformly at the contour of integration. This is satisfied if 



w — \ 
1 — p 



w 



< 1. (5.26) 



The latter is true for a sufficiently large circle, \w\ 3> 1, where this 
expression is arbitrarily close to (1 — p). Summing the geometric 
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progression, we obtain the integral representation (j5.24p for n < N. 
When k becomes negative, keeping the contour large ensures that the 
pole w = 1 turns out to be inside the contour together with w = 0. At 
the moment we also have freedom of putting the pole w = p/{p — 1) 
inside or outside the contour. As this is not a pole for positive r, it 
does not affect the convergence, though changes the part of the sum 
with r < 0. As it will be seen from the proof of orthogonality, it 
should be left outside the contour. 

Now we are in a position to define the functions <&^(r). First we 
note that 

K = Span{T"-\...,l} (5.27) 

and hence can be generated by any n polynomials of powers from zero 
to (n — 1). Indeed the function $"(r) defined by ()5.23p is a polynomial 
of power j. To check the orthogonality relation (j5.2ip we substitute 
the integral representations ()5.24l I5.23P into ()5.2ip and interchange the 
order of summation and integrations. Note that choosing the contours 
as To,! = {\w\ = R} Fi = {\v — 1| = l/R} with R large enough we 
can make the ratio arbitrarily close to (1-p) and hence 



1 — p^w — l)/'W 



1 — p{v — l)/v 



< 1 (5.28) 



uniformly in Fo,i x Fi. Therefore the interchange is allowed. After 
the summation we have 

j; <^nr)^]ir) (5.29) 



tGZ>, 



i+1 



dv f dw \ w ) l^^-i^ rw\^ fwq+p 
J 27ri Ttq ^ 27ri vw{w — v) \v) \vq + p 

where contour Fo,t, contains the pole w = v inside for each v £ Ti, i.e. 
the whole Fi is inside Tq^^. From < k,j < N we see that w = 1,0 
are not poles of the expression under the integral. Therefore we can 
perform the integration in w taking into account the contribution of 
the remaining simple pole w = v, which yields 





(I5.29P = 




,-1^-^- 
V J 


1 

(5.30) 


Making the substitution z 


:={v 


- l)/v, 


we obtain 


the desired orthog- 


onality (I5.2ip 














(|5.29P 


-1 


dz 1, 
z 
27ri 


'-^ = 5..- 


(5.31) 
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Now we note that assumption (A) is fulfilled, 

$^(T)=p = </.n(T„\i,r), (5.32) 

and we can write the kernel. Observe that the function (pnix,y) can 
be written in the form 

(t>n{x,y) = ^A dw—————. (5.33) 



27ri /pjj _^ w[w + 1)^ y 

Several convolutions with itself result in 

dw 



(/)("^'"2)(ri,r2) = l(n2 > ni)p"2-"i 



w 



To.-i 2vri 



(u; + l)^2-ri-i^ (534) 



After a variable change 

w = [ L — p 
we arrive at 



Z-V-' 



1(712 > ni)p (p 



dz 



2TTiz 



2 



ni-n2 / z — 1\'^^^'^'^^^ 

1-p J .(5.35) 

Now let us perform the summation in (j5.22p using the integral 
representations ()5.24l5.23p for ^^(r) and <^^(t) respectively. As the 
integral representation of ^^{t) is identically zero for A; < 0, the sum- 
mation can formally be extended to /c = oo. Then, we bring the 
summation under the integrals, which is allowed provided that 



w^v — 1) 



v{w — 1) 



< 1 (5.36) 



uniformly in Fq,! x Fi. This is satisfied, for example, for Fq,! = {|t(;| = 
R} and Fi = {|u — 1| = 1/R} for R large enough. The contours can 
be continuously deformed without crossing the poles. In particular, 
contour Fi is constrained to be inside contour Fi q. 
Then, the summation yields 

oo 

E^nl-.(n)<^::^.,(r2) (5.37) 

fc=l 

^^ dv J- dw {i-p{^)r{^r{w/vr-^'' 



r, 2TTiv /r„_„ 27riw (i -p(i^))-2+i(i^)n2(y, _ ^) 
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After substituting (|5.35p and (|5.37p into (|5.22p we arrive at the final 
expression for the correlation kernel. The result can be summarized 
as follows. 

Proposition 5.4 The correlation kernel of the measure M, |5.^[j, is 



K{ni,Ti;n2,T2) = (5.38) 

dv r dw {l-p{^)Y'{^T'{w/vY+^ 



P 



r, 27TivJr,^^ 27riw (i -p(-^))-2+i(i^)n.(^ _ ^) 



6 TASEP current correlations 

In present section we obtain the current correlation function of the 
TASEP with step initial conditions and study its asymptotical be- 
haviour. Note that the dynamics of particles is independent of the 
particles on the left. Therefore, one can limit the consideration to 
A^ rightmost particles. The previous results furnish us with all the 
means required to write the quantity of interest ()5.2p as a Fredholm 
determinant, which is the content of the following theorem: 

Theorem 6.1 Given N particles starting at the same moment of time 
t = at sites x^ = 1 — i,i = 1,2, ..., N, let t„4 < ira2 ^ " " " ^ ^n^ 
be the moments of time, when m < N particles labeled by indices 
Hi < n2 < ■ ■ ■ < Um jump off positions Xm, ■ ■ ■ Xn^, respectively, 
where Xi = x + N — i. Then, their joint distribution P, defined in 
h5.S\) . can be represented as Fredholm determinant 



P = det (1 - Xai^Xa))p(|„,,...,„„}xZ>,) ' (6-1) 

where Xa{ni){t) = l(t > a.;). 

Proof. Having written the correlation kernel for the measure M, we 
can use the correspondence stated in the Proposition 15.11 to write 

P = Pro6({t„, <ai}f^{t„2 <a2}p|---p|{t„„ <a„}^ 

= A^(rDK^<ai}f|...f|K™<a„}), (6.2) 

By the exclusion-inclusion principle it is written as a sum 

P = EH^ E E ••• E det{i^(n.„r,;n,^.,r,)}^,,.=i 

?1>0 il,...,in = l ''"l>l^ii Tn>a,i„ 

(6.3) 
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which yields the Fredholm determinant. The formula holds, provided 
that XaKxa is a trace class operator. We do not give a proof of trace 
class property here and refer the reader to similar proofs in other 
papers [D S] ■ ■ 

6.1 Scaling behaviour 

We are interested in the behaviour of the Fredholm determinant in a 
scaling limit. Let us look at a bulk particle with number 

n = uL (6.4) 

jumping off the position 

Xn = xL (6.5) 

at time 

tn = ojL (6.6) 

where L is a large parameter. In the limit L — t- oo the random variable 
ijj becomes deterministic, in a way related to the variables 7 and v. 

According to the law of large numbers for stochastic particle dy- 
namics [261 EH] the density of particles p{x^ t) solves the continuity 
equation 

dp{x,t) dj{x,t) 

where j{x,t) is a stationary state current, which for the TASEP with 
backward sequential update is given by 

i = ?!^^i^. (6.8) 

l-pp 

Solving (j6.7p . (|6.8p with the initial condition p{x,0) = l{—x) we get 
the macroscopic density profile 



pix,t) 



{ 1, x/t<p/(p-l) 

K^-Tt^^)' vKv-\)<xlt<v (6.9) 

0, xji > p. 



At a given time step the number n and the coordinate Xn of a 
particle are related via the density 

I" 00 
n= p{x,t)dx. (6.10) 

Under the scaling ()6.4p - ()6.6p we arrive at the relation 

u = -{^+Vxn^)\ (6.11) 
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which holds when p/{p — 1) < x/^ ^ P- For a rigorous proof of this 
relation one can use [.25j and note that in our case the final positions 
of particles are fixed as Xn = x + N — n. Hence the ratio 

x + N , , 

7=—^, (6.12) 

is fixed as L goes to infinity and therefore 

l = X + v- (6-13) 

This gives the function [25] 

oo{u) = -{^ + ^)\ (6.14) 

where i^ > 0, 7 > 0, which describes the most probable time the n-th 
particle jump off the site x„. 

Under the scaling limit of correlation functions we understand the 
correlations between the jump off times for particles with the numbers 
close to a certain point on the deterministic scale, being apart from 
each other on the scale characteristic for the KPZ class where non- 
trivial scaling behaviour is expected. One chooses |nj — nj\ ~ L^'^. 
Then the dominating scale of time fluctuations is 6ti ~ L^'^. This 
scale defines the domain, where main contribution into the sum ()6.3p 
comes from. 

A rigorous analysis consists of several steps. One has to prove 
that the kernel converges uniformly on bounded sets to its scaling 
limit and that the part of the sum ()6.3p coming from the complement 
to these sets is negligible. We do not follow this program here and 
limit ourselves to a simple saddle point analysis, outlining a sketch 
of the proof of the kernel convergence to the Airy2 kernel. Rigorous 
mathematical details for similar cases can be found in [21[T]. 

Lemma 6.2 Let us introduce scaling variables such that 

m = [L + UiL^/S] (6.15) 

n = [Ll^(1 + niL-i/3) + L^/^Si] (6.16) 

with Ui, Si, i = 1,2 being fixed as L ^ 00, and let K = UKU^^ with 
kernel K given by Proposition \5.3^ and 

U^j{si, s,) = e-W-.("'')+^'^'^'3("''))5ij5(s, - s,). (6.17) 

Then 

lim L^/^k{ni,Ti;n2,T2) = KtKAirySl<.hUl,KtSi;K.hU2,KtS2), (6.18) 

L— >CXD 
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where in the r.h.s. we have the extended Airy kernel, 

KAiry, (6,Ci;6,C2) (6.19) 

^ f /J" dAe^fe-«i)Ai(A + Ci)Ai(A + C2), 6 < 6 

I -/°^(iAe^fe-«i)Ai(A + Ci)Ai(A + C2), 6>6 ' 

and nonuniversal scaling constants 

^' ^ 2(1 + ^)^3(^+^)1/3 ^^-20) 

^(^-1/6^1/6 

''^ ^ (l + V97)2/3(V7 + V^)2/3' ^^'^^^ 



Proof. Suppose 



n, = 


= [Nn] 


Ti = 


= [Aftj(z^i) + Afi/3^ 


a; = 


= [N(l-l)] 



(6.22) 
(6.23) 
(6.24) 



where 



i/i = i/ + niiV-i/3, (6.25) 

and Sj, lij G R, i = 1, 2. As the numbers of particles are constrained to 
be below iV, formally the value of v should be in the range < j^ < 1. 
However, as it will be clear below, this constraint is superficial, and 
can be omitted. 

We also introduce the following functions 

fu{w) = uj{v)\n{q + p/w) + i'\ii{l — \/w) + ^\ii{w) (6.26) 
g{w) = \Ti.{q+p/w) (6.27) 

h{w) = ln(l - l/w) (6.28) 

The position of the double critical point of the function fuiw), which 
satisfies fUwo) = /"{wo) = is 

woiiy) = l + .l^. (6.29) 

V 97 

To obtain the asymptotics of the double integral part of the kernel, we 
recall that it is obtained as a sum J2T=i ^n^-fc(''"i)*^n^-fc(''"2)- There- 
fore we first evaluate the integrals for ^^^_ji.(ri) and ^^^_j^{t2) asymp- 
totically as A^ —7- 00 and then perform the summation. In terms of 
above notations the integrals read 

*S-.(n)) = /J^ („,-i)(,.^rt ("-SD 
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where r = kN^^'^ . 

Instead of the functions in the exponents we use their Taylor ex- 
pansion at the double critical points Wi = WQ{ui), with i = 1, 2 where 
the main contribution to the integral comes from for large N . To 
lowest non-vanishing order we have 

U^{w) = fuM) + \fu{wG){w-w,f + ... (6.32) 

g(w) = g{wi) + g' {wq){w - Wi) + . . . (6.33) 

h{w) = h{w,i) + h' {wq){w - Wi) + . . . (6.34) 

where in the coefficients of the w-dependent terms we, without loss of 
accuracy, replace Vi and Wi by u and wq = wq^v) respectively. We sub- 
stitute these expansion into the integrals, and choose steepest descent 
contours such that they approach the horizontal axis at the points wi 
and W2 at the angles ±7r/3 and ±2-^/3 respectively. Changing the 
integration variables to Ci = {"^ ~ Wi)N^''^f"'{wo)/2 we arrive at the 
integrals defining the Airy functions: 



ooe 



dx ( x^ 



Ai(a:) = / — ^ exp — - xz . (6.35) 

As a result we have the asymptotic behaviour for large N 

f rh'jwp) -sig'{wo) \ 

V (/;(t«o)/2)V3 ; 

^"'-'^^^'^^ ^ {wo-l){qwo+p){NC{wo)/2)y^ ^^-^^^ 
^. f rh'jwp) -s2g'iwo) \ 
^\ (/;(-o)/2)V3 )■ 

Finally, the summation over k can be replaced by an integration 
over r. This requires arguments that the convergence of ^"^„^(ti) 
and $^^_^(t2) is uniform on the sets k < AN^'^ and terms of the sum 
for higher k are negligible for A large enough. These arguments are 
based on the superexponential decay of Airy functions. To perform 
the summations we use one more expansion: 

h{w,i) = h{wo) - h'{wo)wp{u)uiN~^/^ + 0{N~^/^). (6.38) 
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This yields the large A^ behaviour of the sum as 

oo 

Y.KU{ri)KU{r2) = (6.39) 

k=l 

/•oo _ 

X / (iAe^'^'^("2-«i)Ai(A + KiSi)Ai(A + KtS2), 
where 

""^ 2V3 2(07+V79)^/=^(V7 + V^)^/^^ ■ ^ 

^ ^ 2V3ff^(^0) _ p^l/6g-l/6^1/6 

Let us now evaluate the second part of the kernel given by the 
single integral, which can be written as 

/ 27riz^ 
The critical point of the exponentiated function is found to be 

Zc = wq = wq{v) (6.43) 

We use the Taylor expansion around this point to catch the main 
contribution to the integral. Keeping all terms up to the relevant 
order we obtain 

f,^{z)-U,{z) = fuAm) - fu2{w2) (6.44) 



+ 



6N 



f"'{wo){u2-Ui)w'Q{u) ^_^ _ 2 , 

+ 2iVi/3 ^^ WO) +... 



and 



(si - S2)g{z) = sig{wi) - S2g{w2) (6.45) 

+ {U2S2 - uiSi)N~^/^ g' {wq)w' {v) 
+ {si - S2)g' {wq){z - Wq) + . . . 

Substituting these expansions into the integral and integrating along 
the vertical line crossing the horizontal axis at wq we obtain: 

__QQQ -^0 9 

K, (Uj-Ul) (K>,('"l-"2)~'^t(''l-''2)) ~ ~ I ^ 

p 3 4?;ft%2-ui) KhKt(s2«2-Sl«i) 

. (6.47) 
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One can see that the first hne of this expression exactly coincides 
with the factor before the integral in ()6.39p . Furthermore, its expo- 
nential part does not change the value of the determinants, so that it 
can be omitted. The second part can be rewritten using the formula 
from [15] 

^ ^e-{e-f)V4(r'-r)-{r'-r){5+r)/2+{r'-r)3/12 ^g^^g^ 



Y^47r(r' - r 



/oo 
e-A(r-r')Ai(^ + A)Ai(e' + X)dX, (6.49) 
-oo 

where we should set r = KhUi,T' = KhU2,^ = 'i^tsi,^' = KtS2- Notice 
that there is a common factor 

^1,2 = exp(Ar(/^^(^,) _ f,^{w2)) + N^/\sig{wi) - S2g{w2)) (6.50) 

in (6.38) and (6.45). Multiplication of the kernel by the inverse of 
that factor is equivalent to a similarity transformation UKU~^ of the 
operator K with U given by (|6.17p . 

The similarity transformation does not affect the determinants 
det{i^(tij, Sj; Mj, Sj)}i<ij<n and hence the resulting Fredholm deter- 
minant remains unchanged. In addition, after this multiplication the 
limit of the kernel multiplied by A'^^'^ is well defined. As a result, we 
obtain 

lim N'^/^A^\{{Q.'i8) + (6.45)) = KtKAiryo(fihUl,KtSl]KhU2,KtS2)- 

(6.51) 
Finally, we set the large parameter 

L = vN (6.52) 

which removes the dependence of k/j, k^ on v resulting in k^, kj shown 
in ()6.20p , ()6.2ip . The same result could be obtained by setting i^ = 1 in 
the very beginning. The independence of the result on v follows from 
the fact that the dynamics of a particle is independent on particles 
behind it, which is specific of the TASEP. ■ 

The constant kj has occurred before in [141 125) . The Lemma to- 
gether with the steps mentioned above would result in the following 
theorem describing the universal behaviour of the rescaled process 

Theorem 6.3 The following limit holds in a sense of finite- dimen- 
sional distributions: 

hm -— = KtA2{KhU), (6.53) 

L— >-oo L^/'* 
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where A2 is the Airy2 process characterized by multipoint distribu- 
tions: 

Prob(^2(^il) < Si,...,A2iUm) < Sm) 
= det (1 - XsKAny2Xs))L2^{m,...,n^}xB,) ' (6-54) 

7 Conclusion and perspectives 

In the present work we studied multi-point correlations in the dis- 
crete time TASEP with backward sequential update. Using the GGF, 
we have extended the range of analysis to space-time configurations, 
which were not covered by the previous studies. With the help of the 
GGF and the notion of boundaries one can describe the joint probabil- 
ities of both the arrivals at given sites by given times for some particles 
and the hoppings from given sites at given times for the others. The 
arrangement of points in admissible configurations, dealt with by the 
GGF, is defined by the strictly decreasing order of space coordinates, 
(|4.7|) . and weakly increasing order of time coordinates, (j4.8p . The lat- 
ter is exactly opposite to the definition of space-like ordering accepted 
in [SI m [T]. As such, it extends the set of space-time configurations, 
where the probability measure having a compact determinantal form 
can be used as a starting point for calculation of correlation functions. 
However, the constraint ()4.7p is an additional condition, which makes 
the set of point configurations living on the boundaries narrower than 
the whole set of time-like configurations complementary to the set of 
space-like ones. 

We considered the simplest example of the use of the GGF for 
the calculation of correlation functions. This is the case when the 
boundary consists of consecutive vertical lines in the space-time plane 
and the quantity calculated is the joint probability for selected par- 
ticles to jump from corresponding lines by given time moments. The 
expression for the kernel (|5.38p we derived in Section 5 should be com- 
pared with that obtained in |13j for the dynamics of a tagged particle 
in the TASEP with the step initial conditions. Indeed, in the case 
of strictly ordered times ti < ^2 < • • • < ^m , all endpoints in the 
joint distribution correspond to reference points of a hole moving in 
the opposite direction. For the TASEP with parallel update where the 
particle-hole symmetry holds, the hole can be considered as the tagged 
particle. Therefore the difference in the kernels must be attributed to 
the difference in the updates, parallel and backward sequential. One 
also notices a partial similarity between the auxiliary determinantal 
weights obtained in our case, (j5.4p . and proposed in [13j for the case of 
current correlations in the TASEP with parallel update (see formula 
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(6.7) of [13] ). They have similar structure of products of Schur func- 
tions, except one term which seems to have more compUcated form 
in our case. It was also communicated to us by A. Borodin that the 
auxiliary determinantal measure ()5.4p is a particular case of multi- 
parametric family of measures studied in [7] . One of the questions to 
answer is how these parameters can be incorporated into the particle 
dynamics. 

In full generality, the method based on the GGF can be applied 
to obtain the correlation functions on more general sets. There are 
several ways for generalization. One can obtain an auxiliary determi- 
nantal process associated to an arbitrary boundary B. Then, one can 
try to construct a cascade-like determinantal process similar to that 
in [1] using a series of subsequent boundaries. At simplest, these could 
be the fixed coordinate boundaries studied here: 

TV 

gO) = |J{(x(^) -i + N,ti):tie Z}i (7.1) 

1=1 

with x^^' < x^"^' < . . . , the upper index j counting the number of 
points within the correlation functions. The first method would work 
with general admissible configurations of different particles, while the 
second would allow to relax time ordering constraint ()4.8p and to con- 
sider several space-time points for a single particle. Another general- 
ization of interest is an analysis of the TASEP with different initial 
conditions. Note that the GGF allows varying not only initial space 
positions of particles but also the time moments when the particles 
enter the dynamics. 

The asymptotic analysis of the kernel (|5.38|) we obtained shows 
that current correlation functions converge to that of the Airy2 pro- 
cess. This result agrees with the slow decorrelation [12j observed in 
the models of the KPZ class in 1 -|- 1 dimensions. For the same reason 
we expect that a similar result also holds for general boundaries and 
for the configurations accessible for multi-cascade analysis. The most 
intriguing problem is the correlation along the characteristic lines, 
where they are expected to be different from standard KPZ limiting 
processes. A primary analysis shows that the characteristic lines stay 
beyond the domain satisfying the x-ordering constraint (|4.7p . There- 
fore the technique used here has to be sufficiently modified to advance 
in this direction. 
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A Cyclic structure of the determinants 

Go{x,y) and Gi{x,y) 

Lemma A.l Let a = (ui, . . . ,cr]\f) £ Sn be a permutation, x and y 
he two particle configurations from the range Ii2.1^) and 

N 

J'a = '[lF,^_i{x,-y,^,t). (A.l) 

1=1 

1. Then, ift = 0, T^ = unless a is the identity permutation, 

a = {l,...,N) (A.2) 

and 

x = y. (A.3) 

2. Let t = 1, and {ii, . . . ,ik) be a cycle consisting of k > 1 elements 
of a, such that ai-^ = i2,o'i2 = is, . . . ,ai^. = ii. Then, Jv = 
unless up to a cyclic shift of indices 

i2 = ii + 'i-, ■ ■ ■ ,ik = ii + k - I (A. 4) 

and 

Uii — ^ii ) ^i2 — yi2 — ■^*i ^y ■ ■ ■ y ■^ik — Vik — "^U K -\- L. (^A.Oj 

Proof. We first note that 

Fn(x,t) = (A.6) 

unless 

n<x<t for n<0 (A.7) 

and 

X < t ior n > 0. (A.8) 
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For Tfj to be nonzero all the factors in the r.h.s. of (|A.ip have 
to be nonzero simultaneously. Then, from ()A.7p . (|A.8|1 the following 
inequalities hold: 

(Ti - i < Xi - y^i < t, for cn-i <0 (A. 9) 

Xi-yai<t, for (Tj - i > 0. (A. 10) 

Moreover, the particle coordinates vary within the domain (|2.2|) . which 
implies 

i-ai < x„^ - Xi, i - ai < y^. - Vi, ioi ctj - z < 0, (A.ll) 
ai-i < Xi- x^., ai-i <yi-ya^, for ctj - i > 0. (A.12) 

Comparing ()A.9p and (jA.lip we obtain for (i > ai) 

yi<Xi, y^^ < x^^, (A. 13) 

while from (|A.10p and (|A.12p for the case {i < ai) we have 

Xi-yi< t, x„^ - y^^ < t. (A. 14) 

Consider a cycle consisting of A; > 1 elements, (ii, . . . ,ik)- Below 
we adopt the convention i^+i = ii for j = 1,. . . ,N. Let us choose 
r, 1 < r < k, such that the number i^ is the smallest of the numbers 
ii, . . . ,ik- It follows that 

(^ir-i = ir < ir-i, (A.15) 

ir < ir+i = ai^. (A. 16) 

Using ()A.9p and (jA.13p the value of Xi^ can be bounded from above 
and below respectively: 

V^r<Xi,<y^^^,+t. (A. 17) 

1. Let f = 0. The inequality (JA.ITP then yields 

Vir — Vir+n 

which cannot be satisfied together with ()A.16p within the do- 
main ()2.2p . Thus we conclude that for t = only trivial cycles 
with k = 1 are possible for J> to be nonzero. The permuta- 
tion containing only trivial cycles is the identical permutation, 
a = {1,...,N). From 

Fo(x,0) = 5,,o (A. 18) 

we obtain the first statement of the Lemma. 
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2. Let t = 1. First we note that the only way to satisfy ()A.17p is 
to set 

Vir = Xi, = 2/v+i + 1. (A.19) 

This means that the particles with indices ir and ir+i are next 
to each other and are rightmost in the cycle (and in the cluster), 
i.e. 

ir = ir+1 - 1- (A.20) 

Then, applying (|A.9p . (jA.lll and ()A.19p to the case i = i^-i, o-j = 
ir, we obtain 

which suggests that all sites between Xi^ and Xj^_^ are occupied, 
i.e. the particles with indices v and ir-i belong to the same 
packed cluster of particles of the configuration x. 
As the cluster in the configuration x spreads to the left of Xr for 
at least two sites, we have 

Xir+l = Vir + l — ^ir ~ ^- (A. 22) 

If ir-i = ir+1, then k = 2 and the cycle is completed. If A; > 2 
we have 

ir+2 = 0-v+i > ir + 1- (A.23) 

since the particles with indices ir and v+i are rightmost in the 
cluster. We then apply (|A.10p and ()A.14p . which yields 

Xr+2 = Xr - 2 = yr+2 = Ur " 2. (A.24) 

This procedure can be iterated again and again until at k-th. step 
we obtain 

ir+k-l=ir-l- (A. 25) 

The cycle is completed. 
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